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It is generally believed that the old quantum theory, as presented by Niels Bohr in 1913, fails 
when applied to few electron systems, such as the H2 molecule. Here we find new solutions within 
the Bohr theory that describe the potential energy curve for the lowest singlet and triplet states of 
H2 about as well as the early wave mechanical treatment of Heitler and London. We also develop a 
new interpolation scheme which substantially improves the agreement with the exact ground state 
potential curve of H2 and provides a good description of more complicated molecules such as LiH, 
L12, BeH and He2. 



The Bohr model [l| for a one-electron atom played a 
major historical role and still offers pedagogical appeal. 
However, when applied to the simple H2 molecule, the 
"old quantum theory" proved unsatisfactory 0, • Here 
we show that a simple extension of the original Bohr 
model describes the potential energy curves E(R) for the 
lowest singlet and triplet states about as well as the first 
wave mechanical treatment due to Heitler and London 

!• 

We find the Bohr model of H2 admits other solutions 
than the symmetric one he considered (pictured in Fig. 
^) These provide a fairly good description of the 
ground state E{R) (curve 2 in Fig. [3J at large as well as 
small intcrnuclear spacing R, in contrast with the result 
of Bohr (curve 1 in Fig. [3J. 

Clearly the Bohr picture of a molecule goes wrong at 
large R. Any realistic model must show the ground state 
potential energy function dissociating to H+H. Sommer- 
feld, in his seminal book 3], provided an apt assessment: 
"We shall now describe a little more fully the model that 
Bohr has suggested for the constitution of the hydrogen 
molecule H2, although, nowadays, we can take only a 
historical interest in it." After some discussion he asks: 
"But is it correct?" To which he answers: "Only a short 
while ago, even while this book was in its first edition, 
we were inclined to accept it". And later he concludes: 
"Thus the true model of the H2 molecule is still unknown. 
It will hardly be as symmetrically built as the model ex- 
hibited in Fig. 22" . His Fig. 22 is the same as our 
symmetric configuration in Fig. [21 

It is somewhat ironic that the Bohr picture of the 
molecule never caught on. As with the Bohr atomic pic- 
ture, it contains valuable insight, and can provide a good 
analytical description of molecular behavior. Sommer- 
feld even sensed that the symmetric configuration was 
suspect. In Fig. [31 we present a simple continuation of 
the line of thought that Bohr was following which is in- 
deed asymmetric and provides a good quantitative pic- 
ture of H2 at small and large R. We next outline Bohr's 
insightful picture and our extensions. 



Figure [21 displays the Bohr model for a hydrogen 
molecule, in which two nuclei with charges Z\e\ are sep- 
arated by a fixed distance R (adiabatic approximation) 
and the two electrons move in the space between them. 
The model assumes that the electrons move with con- 
stant speed on circular trajectories of radii p± = p 2 = p- 
The circle centers lie on the molecule axis z at the coor- 
dinates Z\ = ±22 = z. The separation between the elec- 
trons is constant. The net force on each electron consists 
of three contributions: attractive interaction between an 
electron and the two nuclei, the Coulomb repulsion be- 
tween electrons, and the centrifugal force on the electron. 
We proceed by writing the energy function E = T + V, 
where the kinetic energy T = p\/2m J rp\/2m for electrons 
1 and 2 can be obtained from the quantization condition 
that the circumference is equal to the integer number n 
of the electron de Broglie wavelengths 2irp = nh/p, so 
that we have T = p 2 /2m = n 2 Ti 2 /2mp 2 . All distances we 
express in terms the Bohr length ciq = fi 2 /me 2 , where m 
is the electron mass, and take e 2 /ao as a unit of energy. 
The Coulomb potential energy V is given by 

z z z z 1 z 2 

V= + — + 17' 1 
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where r a i (i = 1,2) and ru are the distances of the ith 
electron from nuclei A and B, as shown in Fig. [21 (bot- 
tom), ri2 is the separation between electrons. In cylin- 
drical coordinates the distances are 




here R is the internuclear spacing and <fi is the dihedral 
angle between the planes containing the electrons and the 
internuclear axis. The Bohr model energy for a homonu- 
clear molecule having charge Z is then given by (here we 




FIG. 1: Molecular configurations as sketched by Niels Bohr; from an unpublished manuscript 0, intended as an appendix to 
his 1913 papers. 

discuss the case n = 1) of Eq. (J2J. There are four such configurations for which 

P\ = p2 = P, z\ = ±Z2 = z and = 7r, 0; they are 



E = - ( -J + -g- ) +V(p U p2,Zl,Z2,(f>,R). (2) 
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Possible electron configurations correspond to extrema 
I 



pictured in Fig. |21 (upper panel) . 



In Fig. |31 (lower panel) we plot E(R) for the four Bohr 
model configurations (solid curves), together with "ex- 
act" results (dots) obtained from extensive variational 
wave mechanical calculations for the singlet ground state 
and the lowest triplet state, 3 E+ pj|. In the model, 
the three configurations 1, 2, 3 with the electrons on op- 
posite sides of the internuclear axis (</> = tt) are seen to 
correspond to singlet states, whereas the other solution 
4 with the electrons on the same side ((f) — 0) corre- 
sponds to the triplet state. At small internuclear dis- 
tances, the symmetric configuration 1 originally consid- 
ered by Bohr agrees well with the "exact" ground state 
quantum energy; at larger R, however, this configuration 



climbs far above the ground state and ultimately disso- 
ciates to the doubly ionized limit, 2H + +2e. In contrast, 
the solution for the asymmetric configuration 2 appears 
only for R > 1.20 and in the large R limit dissociates to 
two H atoms. The solution for asymmetric configuration 
3 exists only for R > 1.68 and climbs steeply to dissoci- 
ate to an ion pair, H + +H _ . The asymmetric solution 4 
exists for all R and corresponds throughout to repulsive 
interaction of two H atoms. 

The simplistic Bohr model provides surprisingly accu- 
rate energies for the ground singlet state at large and 
small internuclear distances and for the triplet state 
over the full range of R. Also, the model predicts the 
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FIG. 2: Cylindrical coordinates (top) and electronic distances 
(bottom) in H2 molecule. The nuclei Z are fixed at a distance 
R apart. The two electrons rotate about the internuclear axis 
z with coordinates p\ , Z\ and p2 , Z2 respectively; the dihedral 
angle <f> between the (pi, z±) and (p2, 22) planes remains con- 
stant at either <j> — tt or cj> = 0. The sketch corresponds to 
configuration 2 of Fig. [3J with <j> — tt. 



ground state is bound with an equilibrium separation 
R e w 1.10 and gives the binding energy as Eb ~ 0.100 
a.u.= 2.73 eV. The Heitler-London calculation, obtained 
from a two-term variational function, obtained R e = 1.51 
and E B = 3.14 eV 0, whereas the "exact" results are 
R e = 1.401 and E B = 4.745 eV HE|. For the triplet 
state, as seen in Fig. [3J the Bohr model gives remark- 
ably close agreement with the "exact" potential curve 
and is in fact much better than the Heitler-London re- 
sult (which, e.g., is 30% high at R = 2). One should 
mention that in 1913, Bohr found only the symmetric 
configuration solution, which fails drastically to describe 
the ground state dissociation limit. Although a variety of 
modifications were later considered |2j, [lfj , to our knowl- 
edge the other three solutions of the simplest model have 
never been discussed in the literature. One should cer- 
tainly pay tribute to Bohr's planetary model proposed 
long before the development of quantum mechanics. It 
is somewhat ironic that the Bohr model can be derived 
from quantum mechanics in the limit of large dimensions 

m 

We conclude this first portion of our paper with a quick 
sketch of the way the calculations are carried out in order 
to emphasize how simple the present analysis is, as com- 
pared to the many particle Schrodinger equation. For ex- 
ample, for the configuration 2, with Z\ = —Z2 = z, <f> = tt, 
the extremum equations dE/dz = and dE/dp = read 



Z{R/2 



1 symmetric (Bohr) 

Z!=Z2=0, (j>=7t 



2 asymmetric 
/-!=-/. '. , <j)=7t 



3 asymmetric 
Zj=z 2 , (j)=7t 



4 asymmetric (repulsive) 
Zj=-z 2 , (])=0 
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FIG. 3: Energy E(R) of H2 molecule for four electron con- 
figurations (top) as a function of internuclear distance R cal- 
culated within the Bohr model (solid lines) and the "exact" 
ground J E+ and first excited 3 E+ state energy (dots) ||. Unit 
of energy is 1 a.u.= 27.21 eV, and unit of distance is the Bohr 
radius. Please note a similarity between the symmetric con- 
figuration 1 and Bohr's sketch of H2 molecule in Fig. Q 
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which are seen to be equivalent to Newton's second law 
applied to the motion of each electron. Eq. © specifies 
that the total Coulomb force on the electron along the 
z— axis is equal to zero; Eq. specifies that the pro- 
jection of the Coulomb force toward the molecular axis 
equals the centrifugal force. At any fixed internuclear 
distance R, these algebraic equations determine the con- 
stant values of p and z that describe the electron tra- 
jectories. Substituting these values back into Eq. 
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FIG. 4: Energy E(R) of HeH molecule for the shown electron 
configuration calculated within the Bohr model for m = ni = 
ns = 1, Z]f e — 1.954 (solid line) and the "exact" ground state 
energy (dots). 



yields E(R). Similar force equations describe the other 
extremum configurations. 

The simple Bohr model is also useful in describing 
more complicated diatomic molecules. For N electrons 
the model reduces to finding configurations that deliver 
extrema of the energy 




Interpolated Bohr model 



FIG. 5: Ground state E(R) of H2 molecule calculated within 
the interpotated Bohr model (solid line) and the "exact" en- 
ergy (dots) [jj. Insert shows E(R) with no 1/7? term. Curves 
1 and 2 are obtained based on the quantization relative to 
the molecular axis (small R) and the nearest nuclei (large R) 
respectively. Dashed line is the interpolation between two 
regions. 



1 f 

i?=-V^ + nri,r 2l ..,r N ,4 (5) 

In such a formulation of the model there is no need to 
specify electron trajectories and also incorporate non- 
stationary electron motion. One can obtain the energy 
function J5J from dimensional scaling analysis of the 
Schrodinger equation in large-D limit . This provides 
a link between the old (Bohr-Sommerfeld) and the new 
(Heiscnbcrg-Schrodinger) quantum mechanics. 

Next we discuss the ground state potential curve of 
HeH. To incorporate the Pauli exclusion principle one 
can use a prescription based on the sequential filling of 
the electron levels. In the case of HeH the three electrons 
cannot occupy the same lowest level of HcH ++ . There- 
fore, for the configuration with n% = n% = TI3 = I, the 
true ground state energy corresponds to a saddle point 
rather than to a global minimum. Such a configuration 
is pictured in Fig. 0] (insert). In order to obtain the cor- 
rect dissociation limit we assign the helium nucleus an 
effective charge Z^ c = 1.954. Fig. 0] shows the ground 
state potential curve of HeH in the Bohr model (solid 
curve) and the "exact" result (dots) obtained from ex- 
tensive variational wave mechanical calculations || . The 
Bohr model gives remarkably close agreement with the 
"exact" potential energy curve. 

We have found a simple means to improve significantly 
the Bohr model results for bound electronic states. The 
original model assumes quantization of the electron an- 
gular momentum relative to the molecular axis. As seen 



in Fig. |3 this yields a quite accurate description of the 
H2 ground state E(R) at small R, but becomes less accu- 
rate at larger internuclear separation. An improvement 
emerges from the following observation. At large R each 
electron in H2 feels only the nearest nuclear charge. Ac- 
cordingly, as R — > 00, we have two weakly interacting, 
neutral H atoms. Therefore, at large R quantization of 
the momentum relative to the nearest nuclei, rather than 
to the molecular axis yields a better description of the 
physics. This leads to the following expression for the 
energy of the H2 molecule 

1 ( Ti^ 71/^ \ 

E = - I — + — J + V(r a i : r bl ,r a 2,rb2,ri2,R). (6) 

1 \ r al T b2/ 

For m = 7i2 = 1 and R > 2.8 the expression © has 
a local minimum for the asymmetric configuration 2 of 
Fig. |21 We plot the corresponding E(R) without the 1/R 
term in the insert of Fig. (curve 2). At R < 2.8 the 
local minimum disappears and electrons collapse into the 
opposite nuclei. At small 7? we apply the quantization 
condition relative to the molecular axis which yields the 
curve 1 in Fig. [5J To find E(R) at intermediate sep- 
aration we smoothly connect the two regions by a third 
order polynomial (dashed line) . Addition of the 1 /R term 
yields the final potential curve, plotted in Fig. [SJ This 
simple interpolated Bohr model provides good agreement 
with the "exact" potential curve over the full range of R. 

Next we consider the Li 2 molecule. If we neglect inner 
shell electrons of Li then the Li2 molecule becomes similar 
to an excited state of H2 with n\ — ri2 = n = 2 in Eq. JSJl. 
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FIG. 6: Ground state energy E(R) of Li2 molecule calculated 
within the interpolated Bohr model (solid line) and the "ex- 
act" energy (dots). 

Rescaling coordinates in Eq. ijBJ) as r — > n 2 r, R — > n 2 i? 
yields the energy function 
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FIG. 7: Electron configuration and the ground state energy 
E(R) of LiH molecule as a function of internuclear distance 
R calculated within the interpotated Bohr model (solid line) 
and the "exact" energy (dots). 



E 



1 I 1 



1 



n 2 {2r z al 2) 



V(r a i,r bl ,r a2 ,r b 2,ri2,R) 



1,2 



Hence, the ground state potential curve of Li2 can be 
obtained from the ground state E(R) of H2 using the 
following relation 

E Li2 (R) - E Li2 (00) = \ [£h 2 (n 2 R) - £h 2 (00)] . (8) 

The result is shown in Fig. (solid line). For Li2 the Bohr 
model gives the binding energy Eb = 1.10 eV which is 
very close to the "exact" value of E B = 1.05 eV. 

As an example of application of the extended Bohr 
model to other diatomic molecules, we discuss the ground 
state E{R) of LiH. The Li atom contains three electrons 
two of which fill the inner shell. Only the outer electron 
with the principal quantum number n = 2 is important 
in formation of the molecular bond. Applying a similar 
approach to that used to obtain Fig. we find E(R) for 
LiH as shown in Fig. [7\ (solid line), while dots are the 
"exact" numerical answer. This simple extension of the 
Bohr model provides a good quantitative description of 
the LiH potential curve. In this treatment, the essential 
difference from H2 arises simply because in LiH the n = 2 
electron from Li is much more weakly bound than the n — 
1 electron from H, with the result that for LiH the binding 
energy is twofold less than for H 2 and the equilibrium 
separation roughly twice as large. As seen in Fig. [HJ 
the same procedure also gives a good potential curve for 
BeH, a relatively complex five electron system. 
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FIG. 8: Ground state energy E(R) of BeH molecule calcu- 
lated within the interpotated Bohr model (solid line) and the 
"exact" energy (dots). Insert shows the electron configuration 
at large R; only outer shell Be electrons are displayed. 



Finally we show how our very simple analysis yields 
very accurate potential curve for the He2 molecule. We 
apply the Bohr model with momentum quantization rel- 
ative to the nearest nuclei and assume the electron con- 
figuration as shown in the insert of Fig. Then the 
problem reduces to finding minimum of the following en- 
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FIG. 9: Ground state energy E(R) of He2 molecule calculated 
within the Bohr model (solid line) and the "exact" energy 
(dots). Insert shows the electron configuration. 



ergy function 

1 1 2Z 2Z 2 2 
E= ^j + ^j 1 h 

rf r| r x r 2 r\ + r 2 R + r\ - r 2 



R + 2n R-2r 2 R W 

Minimization of this simple expression at fixed R leads 
the potential energy curve pictured in Fig. El (solid line). 
The curve essentially passes through the "exact" dots 
over the full range of R. 

In conclusion, we find a simple extension of the Bohr 
molecular model which gives a clear physical picture of 
how electrons create chemical bonding. At the same 
time, the description is surprisingly accurate providing 
good potential energy curves for relatively complex many 
body systems. 
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